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Abstract 

A Green's function formalism is used to calculate the spectrum of excitations of two neighboring impurities 
implanted in a semi-infinite ferromagnetic. The equations of motion for the Green's functions are determined in 
the framework of the Ising model in a transverse field and results are given for the effect of the exchange coupling, 
position and orientation of the impurities on the spectra of localized spin wave modes. 

Key words: Heisenberg ferromagnet; Ising model-transverse; Green's function; Spin wave; Impurities modes 
PACS: 75.30.Hx, 73.30.Ds, 75.30.Pd 



1. Introduction 

Solids are known to support a variety of elemen- 
tary excitations, such as phonons, polaritons and 
magnons. In order to model the behavior of these 
excitations, one often makes the assumption that 
the media in which they propagate can be treated 
as a perfect crystal with infinite extension. This 
approximation, however, does not alway give a re- 
liable picture of the system. An example is given 
by low-dimension media, such as ultra-thin films, 
in which the presence of surfaces has been shown 
to strongly influence their spectra of excitations. 
That is also true for media containing impurities 
or defects. These different behaviors occur due to 
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the fact that the presence of interfaces or defects 
in an otherwise ideal medium can modify the mi- 
croscopic interactions in the material. Also, the 
low dimensionality of the medium, together with 
the presence of impurities, acts to break the trans- 
lational symmetry of the system, causing signifi- 
cant modifications on the propagation of the exci- 
tations, and also allowing the existence of localized 
excitation modes. 

In the case of magnetic materials, the proper- 
ties of localized excitations such as surface and im- 
purity modes have been widely studied, theoret- 
ically as well as experimentally [1,2,3,4] and sev- 
eral models have been proposed to elucidate the 
dynamics of these modes. Among these, the trans- 
verse Ising model, in particular, has been shown to 
give a good theoretical description of real materi- 
als with anisotropic exchange (e.g., C0CS3CI5 and 
DyPO,}) and of materials in which the crystal field 
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ground state is a singlet [5] . 

The Ising model in a transverse field has been 
applied to ferromagnets to obtain the spin wave 
(SW) spectrum of semi- infinite media [6,7] and 
films [8,9,10]. The SW frequencies associated with 
the presence of a impurity layer in an otherwise uni- 
form semi-infinite medium has also attracted at- 
tention. In both cases the impurities, while break- 
ing the translational symmetry along the direction 
perpendicular to the surface of the media, were as- 
sumed to be uniformly distributed along the plane 
of the film layers. The results showed the presence 
of localized impurity modes which were found to 
depend on the position of the impurity layer within 
the film, as well as on the strength of the exchange 
coupling between the magnetic impurity sites. A 
different aspect of this problem is to consider the 
effect of the presence of localized impurities in the 
medium. In the present paper we make an exten- 
sion of the theory presented in Ref.[ll], by consid- 
ering the effect of several localized impurities in 
an otherwise pure semi- infinite fcrromagnct, in the 
context of the transverse Ising model. Numerical 
results are obtained by means of a Green's func- 
tion technique. The paper is structured as follows: 
in section II the theoretical method is introduced. 
Numerical results are presented and discussed in 
section III. In section IV the main results are sum- 
marized and conclusions are presented. 



2. Model and Green's function formalism 

We consider a semi-infinite ferromagnet with a 
(001) surface and a simple cubic structure (lat- 
tice constant a). Two nearest- neighbor localized 
impurities spins are taken to be embedded in the 
medium at distance (n — l)a from the surface 
(where integer n > 1). The localized spins are 
described by the transverse Ising model and the 
Hamiltonian of the system in the presence of an 
external field can be written as 

ft = ~2 ^2 J i™ s i s ™ - hiS i w 

where Sf and S* are the x and z components 
of the spin operator S, with S = \ for all sites. 



The first term in the right hand side of Eq. (1) 
contains the contribution due to the exchange in- 
teraction. Throughout this paper we assume that 
the summation runs over nearest neighbor sites. 
The second term on the right hand side of Eq. (1) 
refers to the effect of the transverse magnetic field 
in a given site I (this field is h for host sites in the 
bulk and hs for host sites at the surface). For con- 
venience, we shall re-express the Hamiltonian as 
TL = Ho + Hi, where Ho is the Hamiltonian of a 
pure host ferromagnet, whereas Hi corresponds to 
the perturbation caused by impurities. 

Hi = H d + H 'd' + ftoo' (2) 

where the H d (H >d') terms contain the exchange 
coupling between an impurity at a given site la- 
beled o (o') and its neighboring host sites. These 
terms can be written as 

Hod = -(Jo-J)J2 S d S o - ( h » - h ) S o, (3) 

d 

Ho'd' = -(Jo' -J)J2 S d' S o' - i h °' - W" ( 4 ) 
d> 

H 00 > = -(Ji-J)S z S z ,. (5) 

where the d (d 1 ) index labels the nearest neigh- 
bor host sites. The exchange constant assumes the 
values J = J' {J ' = J") for the interaction be- 
tween the impurities and host sites in the interior 
of the medium and J Q = J' s (J > = J$) when both 
sites are at the surface and J otherwise. Likewise, 
the term H 00 ' in the Eq.(2) expresses the exchange 
interaction between the impurities. The Zeeman 
contribution in Eqs.(3) and (4) describes the effect 
of the transverse magnetic field at the impurities, 
h Q and h > , which assume the values h' s and h'g if 
those are located at the surface of the medium, and 
h' and h" otherwise. 

Since we are considering only nearest neighbor 
spin interactions, we can distinguish two physically 
distinct orientations for the impurities, namely: 
one with the two impurities aligned along the x 
axis, which we will refer to as the X case (Fig. la) 
and another, Z case (Fig. lb) with two impurities 
aligned along the z axis. In order to obtain the ex- 
citation spectra of these systems, we extend the 
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Fig. 1. Representation of the interaction scheme for two 
nearest neighbor impurities in a Ising ferromagnct in a 
transverse field. The impurities (black circles) are exchange 
coupled with each other and with their nearest neighbors 
in the host medium, and may be aligned along (a) the x 
axis or (b) the z axis. 

formalism for a single localized impurity presented 
in Rcf.[ll]. Thus, we start by defining the Green's 
functions ((Sf ; S^)) w , where a and (3 stand for 
the Cartesian components of the spin operators 
and w is a frequency label. In the present paper, 
we use the retarded commutator Green's function 
Gj m (q||,a;) = ((Sf ; S^))^. These functions must 
satisfy the equation of motion [12] 

■±(\Sl,S' m ]) (g) 

+ (([Sf,H];Sf n )) 0J . 

Previous calculations for the pure system showed 
that a second order phase transition should oc- 
cur at a temperature Tc, with tanh (h/2kBTc) = 



^((Sf;S z m )) u 



h/3J, where ks denotes Boltzmann's constant, 
such that for T < Tc the average spin orientation 
at each site can have components in the x and y 
directions, whereas for T > T c it lies along the x 
direction. The presence of impurities in the sys- 
tem is expected to change the critical temperature 
to T c , which in some cases may be greater than 
Tc- In order to simplify the calculations, in this 
work we focus on the high-temperature regime 
(i.e. T > Tq). The Green's function for the pure 
system can be obtained by solving Eq. (6) with 
Ti replaced by Tig. The solution, which describes 
an ideal semi-infinite Ising ferromagnet (i.e., with 
translational symmetry parallel to the surface), is 
well known and is given by [13] 

gLM = ~ E G \ M cx p h\ ' ( r; - r ")]' 

q 

(7) 

where the vectors rj and r m indicate the positions 
of two given sites I and m, q|| = (q x , q y ) is an in- 
plane wave vector and M is the number of sites 
in any layer parallel to the surface. The Fourier 
amplitudes of this function are found as 

G «" {u>)= 2irJ( X 1 -x-i) 

x | x 1 "-™' 1 - 1 + X A A x"+"' 
1 + xA 



(8) 



Here the labels n and n' are layer indices for the 
sites / and m, respectively (with n = 1 being the 
surface layer.), and x is a complex number that 
satisfies the condition 

x + x- 1 = [h 2 - Ah JR x i(c\\\ ) - uj 2 ]/hJR x (9) 

where 7(qii) is a structure factor that, in the case 
of a simple cubic lattice, is given by 7(qy) = 
T;{cos(q x a) + cos(q y a)]. In addition, when describ- 
ing localized modes the condition \x\ < 1 must 
also be fulfilled. The parameter A in Eq. (8) con- 
tains the information regarding the surface and is 
given by 



A = 



Lo 2 (h s R x s -hR x )- hh s (hR x - h s R x ) 



hh s R x R x s J 



4 7(qn) 



ii. 
J 



(10) 
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where, from mean field theory, we obtain the spin 
averages at the bulk and surface as R x = (S x ) = 
tan\\{h/2k B T) and R% = tanh{h s /2k B T) 1 re- 
spectively. 

The presence of localized impurities in an other- 
wise ideal medium acts to break the translational 
invariance of the system. Consequently, the calcu- 
lations for the impure system must be performed 
in real space. By including the effects of the impu- 
rity in the Hamiltonian and applying Eq. (8), one 
can obtain a new Green's function Gi m (u>), which 
is found to obey the equations 

AijGi m (u) = Sij- 
2ir 

— [P^-U^-U'^-U'^GUu) 

(11) 

where 

_ uj 2 -[h+ (h - h)5 lo + {h , - h)8 lo ,] 2 
lj ~ h+[h+(h -h)5 l0 + (h , -h)6 w ] u ' 

p 

Ulj =^2(Jo- J)R?[$lo5dm+5ldSjo}, 

d 

U 'ij = E( J °' - J)RilSio'6 jd ' + S ld >6 jo ,}, 

d> 

C//j = (J/-J) R![5io>6 jo + 6io6 jo ,}. 
where the index o (o') assumes the values (2) for 
each impurity site in the X case and (6) in the 
Z case (see Fig. 1). The second summation runs 
over 5 neighboring host sites to either impurity. By 
rewriting Eq. (11) in matrix form we obtain the 
Dyson equation 

[(G°H)- 1 -V]GH = I, (12) 

where G°(u>) and G(w) are square matrices 
with elements given by (2ir/Rf )G° m (w) and 
(2ir/Rf)Gi m (Ld), respectively. I is the unit matrix, 
V is an effective potential related to the impurity 
term TLi, with elements 

_L0 2 -h 2 

Vlj ~—f^% - A u - (13) 

-Uu-U'v-U'Sy 

Thus, Eq. (9) can be said to relate the matrix 
Green's function of the pure (G°(w)) and impure 
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Fig. 2. Localized spin wave frequencies as a function 
of the exchange parameter J , for the X case, with 
h a = h a i = 0.65/i and T = 2.5J/kg. The parameters are 
J i = 1.5 J and Jj = 0.25J when both impurities are at the 
surface (solid lines), and J a i = 2.5J, Jj = 0.25J when the 
two impurities are at the second(dashed lines) and third 
layers (dotted lines). 

(G(cij)) systems. The spectrum of localized modes 
is then found by numerically calculating the fre- 
quencies that satisfy the determinantal condition 

det[I-G°(w)V] =0 (14) 

which gives the poles of the Green's function for 
the impure system. 



3. Numerical Results 

Impurities modes spectra were calculated as 
functions of the exchange and effective field pa- 
rameters. In order to assess the influence of the 
position of each impurity on the excitation spec- 
tra, we obtained numerical solutions of Eq.(14) for 
impurities located at the surface of the medium 
(layer 1) as well as deeper in the system (layers 2 
and 3). Specifically, for the X case, we considered 
three configurations, namely, we set n = 1 (both 
impurities at the surface), n = 2 (both impurities 
at the second layer) and n = 3 (both impurities at 
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Fig. 3. Localized spin wave frequencies as a function of 
the exchange parameter J in the upper impurity in an 
Ising ferromagnet in a transverse field, for the Z case, 
with h = h a i = 0.65fe and T = 2.5J/kg. The parameters 
are J i = 1.5J and Jj = 0.25J when the upper impurity 
is located at the surface (solid lines), and J a i = 2.5J, 
Jl = 0.25J when the upper impurity is at the second 
(dashed lines) and third layers (dotted lines). 

the third layer) according to the notation in Fig. 
la. Likewise, for the Z case, following Fig. lb we 
set n = 1 (i.e. one of the impurities at the surface, 
the other at layer 2), n — 2 (one at layer 2, the 
other at layer 3) and n = 3 (impurities at layers 3 
and 4). 

Figure 2 shows the local SW impurities spectra 
for the X function of the exchange con- 

stant for the coupling between one of the impuri- 
ties and its neighboring host sites (J D ), while the 
parameters for the second impurity are kept con- 
stant. We have used J > = 1.5 J and Jj = 0.25 J for 
n = 1, and J > = 2. 5 J, Ji — 0.25 J for n = 2 and 
3. The field parameters were h a = h > = 0.65h and 
we used T = 2.5J/fc B in all cases. The three sets 
of branches correspond to the impurities located in 
layers 1 (solid line), 2 (dashed line) and 3 (dotted 
line). The impurities modes can be classified as res- 
onance modes (i.e. those occurring in the SW bulk 
band) and defect modes (those found outside the 
bulk band). In this work we consider defect modes, 




Jo,/ J 

Fig. 4. Localized spin wave frequencies as a function of 
the exchange parameter J a i in the lower impurity in an 
Ising ferromagnet in a transverse field, for the Z case, with 
h = h D i = 0.65ft and T = 2.5J/ks. The parameters are 
Jo = 1.5J and Jj = 0.25J when the upper impurity is lo- 
cated at the surface (solid lines), and Jo = 2.5J, Ji = 0.25J 
when the upper impurity is at the second (dashed lines) 
and third layers (dotted lines). 

since these are easier to measure. Therefore we re- 
strict the calculation to the region outside the bulk 
SW band. 

The lower limit of the bulk band is represented 
in the graph by the horizontal line. In the absence 
of the coupling between the impurities, the graph 
would show a horizontal line (for the impurity with 
J' constant) superimposed on the decaying curve 
of the SW branch associated with the second im- 
purity, which merges with the bulk band for low 
values of J Q . The exchange coupling thus creates 
a mode repulsion effect when J D w J >. As ex- 
pected, the magnitude of this effect depends on the 
strength of the exchange coupling between the im- 
purities. When both impurities are at the surface, 
the lower coordination number causes a large fre- 
quency shift in comparison with the other config- 
urations. 

The effect of the alignment of the impurities is 
shown in Fig. 3, which presents results for the SW 
frequencies in the Z a function of the ex- 
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change constant at the upper impurity, while the 
remaining parameters were kept constant. For n = 
2 and n = 3, the results are similar to the ones 
obtained for X case, apart from a small frequency 
shift. For n = 1, however, the branches are strongly 
shifted. These results demonstrate that the influ- 
ence of orientation becomes particularly important 
for impurities close to the surface. This is a con- 
sequence of the smaller number of neighbors, to- 
gether with modified exchange parameters at the 
surface. 

Figure 4 shows results of frequency as a function 
of J i at the lower impurity in the Z case, with 
the remaining parameters kept constant. In this 
case, the exchange parameter being varied is asso- 
ciated with the lower impurity (see Fig. 1), which 
was assumed to be located at the second, third and 
fourth layers. In contrast with the results of Fig. 
3, the frequency branches are not found to merge 
with the bulk band for low values of the exchange 
parameter. Moreover, when the lower impurity is 
localized at the second layer (solid line), the curves 
show a larger shift in comparison with the remain- 
ing configurations. This effect is a consequence of 
the modification of the exchange parameters for 
the upper impurity when located at the surface. 

The influence of the local fields is shown next. 
Figure 5 shows the local SW frequencies as a func- 
tion of local field parameter h at one of the impu- 
rities, for the X case. The remaining parameters 
are the same as in the previous graphs. The results 
were obtained for n = 1 (solid lines) ,n = 2 (dashed 
lines) and n = 3 (dotted lines). The graphs show 
the mode repulsion, due to the coupling of the im- 
purities. Also evident is a large frequency shift for 
n = 1, due to the modified exchange parameters 
at the surface. For fields larger than h a = 0.69 h 
(n = 1), h a = 0.99 h (n = 2) and h Q = 1.11 ft 
(n = 3), the upper frequency branches are found 
to merge with the bulk band, thus becoming reso- 
nance modes. 

A similar graph is shown in Fig. 6, this time 
for the Z case. The results were calculated for a 
varying local field at the upper impurity. In con- 
trast with the previous figure, the higher frequency 
branches in the graph do not display any notice- 
able frequency shift for low values of the local field. 
On the other hand, the high frequency excitations 
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Fig. 5. Localized spin wave frequencies as a function of 
the local field parameter h in one of the impurities, in 
an Ising ferromagnet in a transverse field, for the X case. 
The remaining parameters are the same as in the previous 
graphs. The results were obtained for both impurities at 
the first (solid lines), second (dashed lines) and third layers 
(dotted lines). 

are observed to become resonance modes for val- 
ues of local fields that approximately the same as 
in the X case. 

Figure 7 shows a plot of frequency as a function 
of the exchange coupling constant for the interac- 
tion between the impurities (J/). The graph shows 
results for both the X and Z cases, and the re- 
sults were obtained for ferromagnetic ( Ji > 0) as 
well as antiferromagnetic (J/ < 0) coupling. For 
the X case, we assumed that both impurities were 
located at the surface of the system, whereas in 
the Z case we considered the upper impurity at the 
surface, and the lower one at the second layer. The 
branches behave quite distinctly in each case, with 
the results for the Z case displaying a large shift 
in comparison with the other case. Also prominent 
in the X case is a mode crossing that is found for 
Ji = — 0.86J. In contrast, the results for the Z 
case show no mode crossings, although the low fre- 
quency is observed to have a maximum at Jj = 
— 1 . 1 U and the high frequency branch has a min- 
imum at Ji = —1.37 J. 
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Fig. 6. Localized spin wave frequencies as a function of the 
local field parameter h for the upper impurity, in an Ising 
fcrromagnct in a transverse field, for the Z case, when the 
upper impurity is located at the first (solid lines), second 
(dashed lines) and third layers (dotted lines). 



4. Conclusions 

We have presented a Green's functions calcu- 
lation of the SW frequencies of localized modes 
associated with two coupled localized impurities 
implanted in an otherwise ideal ferromagnet. The 
results were obtained in the context of the Ising 
Model in a transverse field and show the influ- 
ence of the exchange coupling of the impurities on 
the excitation spectra of the system. This coupling 
modifies the spectra in relation to the results for 
single impurities, especially when the coupling pa- 
rameters for the interactions between the impu- 
rities and the host sites are of the same magni- 
tude. The results also point to a strong effect of the 
orientation of the impurities on the localized SW 
frequencies, especially when the impurities are lo- 
cated close or at the surface. Further studies could 
investigate the effect of the presence of the impu- 
rities on thin ferromagnetic films, as well as the ef- 
fect of their coupling and orientation on the critical 




Fig. 7. Localized spin wave frequencies as a function of 
the exchange parameter J/ for the coupling between the 
impurities. The graph shows results for both the X case 
(solid line) where the impurities are located at the surface 
layer, and the Z case (dashed line) where one impurity is 
located at the surface and the other in the layer below. 
The remaining parameters are the same as in the previous 
figures. 

properties of the system as well as on the overall 
magnetization of the material. 
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